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Recent experiments on polarized deep inelastic lepton-nucleon scattering from SMC 
and E143 have confirmed the finding of the earlier EMC results that the flavor- 
singlet is small 0. In the deep inelastic limit, the integral of the polarized structure 
function is related to the forward matrix elements of axial currents from the operator product 
expansion Combined with the neutron and hyperon decays, the flavor- singlet axial 
coupling is extracted. Since the axial current is the canonical spin operator, g'^ is thus 
the quark spin content of the nucleon; i.e. g^^ = Au + Ad + As, where the spin content 
Aq{q = u, d, s) is defined in the forward matrix element of the axial current, (ps|gi7^75g|ps) = 
2MNSf,Aq. 

The fact that g% which represents the quark spin contribution to the proton spin, is 
found to be much smaller than the expected value of unity from the non-relativistic quark 
model or 0.75 from the SU(6) relation ( 3/5 of the isovector coupling g\ = 1.2574) has 
surprised physicists who have taken the valence quark model for granted. In view of the 
exciting puzzle presented by the experiment, the gross failure of the valence quark model, 
and its possible manifestation of the axial anomaly besides vr*^ —>■ 77, several lattice QCD 
calculations were carried out in the last 6 years as an attempt to meet the challenge of 
explicating the phenomenon directly from QCD. 

In the following, I shall summarize the progress that has been made in this direction and 
point out the pitfalls and/or shortcomings in each of the calculations reviewed here. There 
are three approaches to address the issue of the proton spin components. The first is the 
calculation of the gluon spin content, the second is the calculation of the axial anomaly, and 
the third is the calculation of the quark spin content from the axial- vector current. 

1 Gluon Spin Content 

Ref 0: Jeffrey E. Mandula, "Lattice Simulation of the Anomalous Gluon Contribution to 
the Proton Spin", Phys. Rev. Lett. 65, 1403 (1990). 

This is the first attempt to calculate the gluon spin content as a way of understanding 
the EMC experiment. 

1.1 Calculational details 

• Lattice parameters: 204 6^x 10 lattice at (3 = 5.7 and the Wilson mass parameter 
K = 0.162. 

• This a quenched approximation (without dynamical fermion loops) with lattice spacing 
a ~ 1.0GeV-\ 

• The gluon spin operator is the lattice version of the anomalous current. 

= e^^xpTrA^^G^P - 2/?>A^AP) (1) 
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which is a particular case of the Chern-Simons form [|^. The gluon spin fraction is 
defined as 

{ps\K^\ps) = 2MNS^Ag. (2) 

Since the current is not gauge invariant, the calculation is done in the ^4 = gauge 
where the second term in Ki (i = 1,2,3) vanishes. 

It turns out that the signal is very noisy and only an upper bound is given, i.e. Ag < 
1 0.5 1 . Since the original motivation of the work is to evaluate the possible gluon spin 



mixture in the quark axial-current matrix element via the triangle diagram g, |T0 
the author gives the upper limit of the mixture {3as/27i)\Ag\ < 0.05. 



1.2 Comments 

• This work has drawn several comments : 

Efremov, Soffer, and Tornqvist Jll]! assert that what is calculated in ^ contains a 



contribution from the "ghost" pole. On the other hand, Manohar [|T3| points out that 
even in the forward direction, Ag may still depend on the gauge direction rj. Mandula 
replies [^] that "while gauge-variant singularities and rj dependence are kinematically 
allowed in the forward matrix element {ps\Kf^\ps), they are not expected to occur, on 
either dynamical, symmetry, or invariance grounds." 

We believe that the issue of "ghost" pole is a red herring. The "ghost" was invented 
by Veneziano []T^ for the convenience of taking care of a minus sign in the topological 



susceptibility which may well come from a contact term. Unlike the claim of Efremov, 
Soffer, and Tornqvist that "This ghost pole is necessary for the resolution of the U(l) 
problem in QCD", Witten's resolution of the U(l) problem in terms of the topological 
susceptibility does not depend on the "ghost". Furthermore, it has been shown 



that one can adapt Veneziano's approach without invoking the "ghost" . On the other 
hand, Manohar does have a point, especially in view of the fact that it has been shown 
by Cronstrom and Mickelsson |^ that under gauge transformations. 



The last term is not sensitive to the local deformation of U and can be considered a 
topological current. For large gauge transformation where the topological (or instan- 
ton) number changes, the forward matrix element 

\im{p + q\d''H^,{U)\p) = \iimq''{p + q\H^,{U)\p) (4) 



does not vanish |]T5|, |T^. Thus the calculation of the gluon spin content using the 



anomalous current in eq. (|l]) is gauge-dependent even for forward matrix element. 

It turns out that in the temporal axial gauge, ^4 = 0, the anomalous current Ki 
in eq. (0) coincides with the gluon spin operator 2Tr{A x Ey Hence, modulo 

gauge-dependence, what is calculated in is the gluon spin content. 
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2 Flavor- Singlet g\ from Axial Anomaly 

In this approach, one calculates the following matrix element 




The basis for calculating A instead of the matrix element of the axial current directly 
can be seen as follows. Inserting the anomalous Ward identity 



Nf 

d^Al = 2 Y: rUfq^H^qf + Nj-^Gl^G^^^ (6) 
/=i 

between proton states, one obtains 



u{p)il^u{p')[2Mg\{q^) + q^h\{q^)] =< p\2mP + 2Njq\p' >, (7) 

where P is the pseudoscalar current, i.e. mP = YJfLi^f^fil^^lf and q is the topological 
charge operator, q = (as/47r)G'^^G"''^''. 

Since there is no Goldstone pole at the chiral limit in the U(l) channel, one can take the 
— * and the m —>■ limits and drops the and P terms in eq. (0) to obtain 



= AS = A (8) 

Thus one can calculate the matrix element of the axial anomaly to obtain AS. One will 
not be able to calculate Au, Ad, and As separately in this approach though. 

2.1 Quenched Approximation 
2.1.1 Two calculations 



1. Ref [|1^]: B. AUes, M. Campostrini, L. Del Debbio, A. Di Giacomo, H. Panagoulos, 
and E. Vicari, "The Proton Matrix Element of the Topological Charge in Quenched 
QCD", Phys. Lett. B336, 248 (1994). 

• Lattice parameters: 100 16^x32 lattice at j3 = 6.0. Smeared source for the nucleon 
in the Coulomb gauge. Wilson hopping parameters k = 0.153,0.154,0.155 are 
used for the quark propagators. 

• Nonperturbative method is used to determine the finite lattice renormalization 
constant. 

• The signal is zero within errors. An upper bound of the matrix element A in eq. 
(§) is set to be 0.4. 
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2. Ref |]20|: R. Gupta and J.E. Mandula, "Matrix Elements of the Singlet Axial Vector 
Current in the Proton", Phys. Rev. D50, 6931 (1994). 



• Lattice parameters: 35 16^ x 40 lattice aX (3 = 6.0. The quark masses used 
are k = 0.154 and 0.155, corresponding to pions of mass 700 and 560 MeV, 
respectively. The Wuppertal source is used for the nucleon, hence it is a gauge 
invariant calculation. 

• Again, the authors found that the data are statistically indistinguishable from 
zero. 



2.1.2 Comments 

• Irrespective of the fact that no signal is seen in these calculations, it is pointed 
out emphatically by Gupta and Mandula that this approach which relies on 
the use of the axial anomaly fails in the quenched approximation. This can be 
see as follows. What is derived in eq. (H) depends on the fact that there is no 
massless Goldstone mode in the U(l) channel at the chiral limit. However, this 
is not realized in the QUENCHED approximation where the would be Goldstone 
bosons are present. It is clearly seen in the study of the anomalous Ward identity 
for g\ 1 18] that the induced pseudoscalar form factor Ha has a pole structure 



where m^^ is the would-be Goldstone meson mass, while the pseudoscalar 

^ 2 

TTl 

current 2mP has a pole structure ^2_^2 ■ Combined, they give the structure 

'^2_™2° which tends to unity at the — > and m — > limits. In other words, 

in the quenched approximation, g'^ is not simply related to A as in eq. (|). The 
pseudoscalar current 2mP also contributes in this case even in the chiral limit as 
is in the Goldberger-Treiman relation. 



Gupta and Mandula [g^] further suggest that due to the presence of the dou- 
ble/single would-be Goldstone pole in the pseudoscalar current /axial anomaly, 
the result from the quenched approximation diverges at the chiral limit and thus 
one will not be able to utilize the anomalous Ward identity to calculate g^. For- 
tunately, it can be seen from our earlier study []T^ that the double pole in the 
pseudoscalar current cancels the single pole from the axial anomaly. This is due 
to the fact that the topological susceptibility vanishes at the chiral limit so that 
the vacuum has no 6 dependence. Furthermore, the remaining single would-be 
Goldstone pole (at = 0) will be cancelled by the quark mass in the pseu- 
doscalar current 2mP. Thus every term in the expression of eq. (|^) is finite in 
the chiral limit of the quenched approximation. We conclude that one can still 
utilize the anomalous Ward identity in the quenched approximation, provided the 
pseudoscalar current is evaluated in addition to the axial anomaly term. 
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2.2 Calculation with Dynamical Fermions 

Ref. [|2T|: R. Altmeyer, M. Gockler, R. Horsley, E. Laermann, and G. Schierholz, "Axial 
Baryonic Charge and the Spin Content of the Nucleon: a Lattice Investigation" , Phys. Rev. 
D49, R3087 (1994). 

2.2.1 Calculational Details 

• Lattice parameters: 16^ x 24 dynamical fermion lattice at (3 = 5.35 with four flavors 
of staggered fermions with mass ma = 0.01. 85 gauge configurations separated by 5 
trajectories are used. The lattice spacing determined by the p mass is 0.14 fm. 

• The result of eq. (D calculated at |gl ~ 500 MeV gives A = 0.18 ± 0.02 which the 
authors identify with AS through eq. (|). 

2.2.2 Comments 

• This is the first calculation which uses the dynamical fermion gauge configurations and 
actually sees a signal which is quite commendable. Moreover, the result agrees with 
the experimental finding. 

• However, this calculation is done at finite q and finite quark mass m. As we alluded 
to earlier, the induced pseudoscalar form factor Ha will contribute at finite q (NB. 
eq. (1^) and the pseudoscalar current 2mP will contribute at finite m. As long as the 
quark mass (in both the valence and the fermion loops) and the momentum transfer 
q are not extrapolated to 0, the relation in eq. (||) will always be subjected to these 
systematic errors. In this approach, one can not separately calculate Au, Ad, and As. 
Furthermore, just knowing through the matrix element of the anomaly does not 
teach us why g\ is as small as it. 

3 Quark Spin Content from the Axial- Vector Current 

This is the direct way of calculating the quark spin content, since it is defined by the forward 
axial current matrix element and one can calculate it flavor by flavor. The only problem is 
that in addition to the connected insertion, there is a disconnected insertion. The discon- 
nected insertion requires the calculation of quark loops which is prohibitively time-consuming 
if one is to invert the whole quark matrix directly. This is the reason why earlier attempts 
were made to calculate the axial anomaly in order to circumvent this numerical difficulty. It 
turns out that there have been three calculations which tackle this problem with different 
numerical techniques. 
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3.1 Quark Loops in a Smaller Box 



Ref. |22]: Jefferey E. Mandula and Michael C. Ogilvie, "A New Technique for Measuring 



the Strangeness Content of the Proton on the lattice", Phys. Lett. B312, 327 (1993). 



3.1.1 Technical Details 

• Lattice parameters: 16 16^ x 24 quenched lattice at /3 = 5.7. The valence quark masses 
are about 0.6, 0.2, and 0.1 GeV corresponding to Wilson k = 0.140,0.160, and 0.164. 
The strange mass in the loop is about 0.16 GeV (k = 0.162). 

• The quark propagator for the loop is calculated in a 9^ box and the loop is evaluated 
on a time slice for every other site in each spatial direction. This saves a factor ~ 500 
in computer time comparing to the straight forward calculation and thus makes the 
project feasible. 

• As tends to be negative for light valence quarks in the nucleon. However, no plateau 
is reached to give a final result. 



3.1.2 Comments 

This technique should work. The devil is in the statistics. 



3.2 Volume Source Technique 



Ref. ||2J]: M. Fukugita, Y. Kuramashi, M. Okawa, and A. Ukawa, "Proton Spin Structure 
from Lattice QCD", Phys. Rev. Lett. 75, 2092 (1995). 



3.2.1 Technical Details and Results 

• Lattice parameters: 260 16'^ x 20 quenched lattice at (3 = 5.7. The Wilson hopping 
parameters n = 0.160,0.164, and 0.1665 are used. The strange quark corresponds to 
K = 0.160. 

• The quark loops are calculated with the volume source technique which relies on 
the cancellation of gauge non-invariant part through gauge averaging. The time slice 
from which the nucleon source emerges is fixed in Coulomb gauge. 

• The results are summarized in the following table. 



3.2.2 Comments 



• The result of AE agrees with experiments well and its smallness is understood as due 
the sea-quarks which give negative contributions to AE. 
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Table 1: Quark spin contents of the proton form the volume source calculation 



g'iiAE) 0.18(10) 

gl 0.985(25) 

Au 0.638(54) 

Ad - 0.347(46) 

As - 0.109(30) 

Fa 0.382(18) 

Da 0.607(14) 

Fa/Da 0.629(33) 



g\, Fa, and Da are smaller than the experiments by ~ 25 %. Possible sources of 
systematic errors include the scaling violation effect due to a fairly large lattice spacing 
Q and the fitting procedure with the ratio method which can give a jackknife bias of 
^ 10% as compared to the more desirable simultaneous fitting procedure p5[. 



• Since the time slice from which the nucleon source emerges is fixed to the Coulomb 
gauge, the results are gauge dependent. Furthermore, the quark current operator 
can mix with the gauge non-invariant anomalous current K^^ in eq. ([l|), a problem 
deja vu. To see this, one can construct the lattice version of i^T^ with staples in such a 
way that their ends touch the time slice which is gauge fixed or use the version given 
by Mandula |Q with the sandwiching U-links lie in the gauge-fixed time slice. These 
Kfj^ matrix elements do not vanish on this lattice and become a part of the results 
obtained in this work. It is then necessary to disentangle this gluon spin contribution 
from the quark spin part before addressing the proton spin components. 

3.3 Stochastic Estimation with Zo Noise 



Ref. p6|: Shao-Jing Dong, Jean- Francois Lagae, and Keh-Fei Liu, "Flavor-Singlet gA from 



Lattice QCD", Phys. Rev. Lett. 75, 2096 (1995). 

3.3.1 Technical Details and Results 

• Lattice parameters: 24/50 16^ x 24 quenched lattice configurations aX (3 = 6.0 for the 
connected/disconnected insertion. The Wilson hopping parameters are k = 0.148, 0.152, 
and 0.154. The strange quark corresponds to k = 0.154. 

• No gauge-fixing is done. Hence the calculation is gauge invariant. The quark loops are 
calculated by the stochastic method with the optimal Z2 noise which is shown to 



give the correct answers for the forward matrix elements of the vector and pseudoscalar 
currents 
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• Finite ma correction for the quark loops of the Wilson fermion is taken into account 

• The results are summarized in the following table. 

Table 2: Axial coupling constants and quark spin contents of proton in comparison with 
experiments 

This Work Experiments 



g'^^ = Au + Ad + As 


0.25(12) 


0.22(10) i / 0.27(10)0 


g\ = Au- Ad 


1.20(10) H 


1.2573(28) 


g% = Au + Ad- 2 As 


0.61(13) 


0.579(25) [gg 


Au 


0.79(11) 


0.80(6)i /0.82(6)i 


Ad 


- 0.42(11) 


- 0.46(6)0 /- 0.44(6) § 


As 


- 0.12(1) 


-0.12(4)i /- 0.10(4) i 


Fa = (Au - As)/2 


0.45(6) 


0.459(8) 


Da = {Au-2Ad + As)/2 


0.75(11) 


0.798(8) 


Fa/Da 


0.60(2) 


0.575(16) H 



3.3.2 Comments 



This calculation is gauge-invariant in that no gauge-fixing is applied. Since there is 
no gauge-invariant dimension-three axial operator for the gluon spin as we discussed 
earlier [14, 1^, these results are not mixed with the gluon spin. Only the quark spin 



content contributes. 

The physical picture of g\ is starting to come in view. The smallness of the quark 
spin content compared to the non-relativistic value of unity is, first of all, due to the 
fact that the combined relativistic effect and polarization of the cloud-quarks reduces 
the connected insertion to 0.62 ± 0.09, a value very close to g\, i.e. g\c,on — 9a- This 
is because the disconnected insertion is almost independent of the flavors u, d, and s. 
Furthermore, the sea-quark polarization is large and in the opposite direction of the 
proton spin. It is the sum of all these effects that produces a small g^. 

It is interesting to observe that the sea-quark spin contribution is independent of the 
quark mass in the loop within errors. This is reminiscent of the 75 current in the 
context of the topological charge and susceptibility. It is suggestive of the importance 
of the zero modes and the instantons. 

Albeit the results all agree with experiments, the systematic errors due to the finite 
volume, discretization, renormalization, and quenched approximation (which could be 
as large as 7% - 20%||3^) need to be addressed before a final precise comparison with 
experiments can be made. 
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3.4 Moments of Polarized Structure Functions 



31 



Ref. 

A. Schiller 



M. Gockler, R. Horseley, E.-M. Ilgenfritz, H. Perlt, P. Rakow, G. Schierholz, and 
"Polarized and Unpolarized Nucleon Structure Functions from Lattice QCD", 



DESY 95-128, |liep-lat/ 9508004] . 



3.4.1 Technical Details 

• Lattice parameters: 400 - 1000 16^ x 32 quenched lattice configurations at (3 = 6.0. 
The hopping parameters are k = 0.1515,0.153, and 0.155. 

• gA = 1-07(9), the connected insertion part of (/^ is 0.59(5). The gi sum rule is 0.166(16) 
for the proton and -0.008(09) for the neutron. 



3.4.2 Comments 

• The results are from the connected insertion only. They have very high statistics. The 
results on the disconnected insertions are forth-coming. 

4 Future 

So far, the lattice calculation is only beginning to shed some light on the quark spin content 
part of the proton spin. We need to calculate the orbital angular momentum and the gluon 
spin in order to complete the picture of the proton spin composition. In the future, the gauge 
invariant calculation of the axial-vector matrix elements to address the infinite volume and 
continuum limits and the dynamical fermions is needed to study the systematic errors of the 
quark spin content. 

This work is partially supported by DOE Grant DE-FG05-84ER40154. 
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